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Abstract 

Low-energy dipole excitations in neon isotopes and = 16 isotones are calculated with a 
fully consistent axially-symmetric-deformed quasiparticle random phase approximation (QRPA) 
approach based on Hartree-Fock-Bogolyubov (HFB) states. The same Gogny DIS effective force 
has been used both in HFB and QRPA calculations. The microscopical structure of these low-lying 
resonances, as well as the behavior of proton and neutron transition densities, are investigated in 
order to determine the isoscalar or isovector nature of the excitations. It is found that the = 16 
isotones ^^O, ^^Ne, ^^Mg, and ^*^Si are characterized by a similar behavior. The occupation of 
the 2si/2 neutron orbit turns out to be crucial, leading to nontrivial transition densities and to 
small but finite collectivity. Some low-lying dipole excitations of ^^Ne and ^"^Ne, characterized by 
transitions involving the z^ld3/2 state, present a more collective behavior and isoscalar transition 
densities. A collective proton low-lying excitation is identified in the ^'^Ne nucleus. 

PACS numbers: 21.60.Jz, 21.10.Gv, 27.30.+t 
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I. INTRODUCTION 



Since the beginning of the nuclear physics, dipole excitations of nuclei represented one 
of the most important ways to investigate nuclear structure. This is due to the fact that 



they can be induced by 7 radiation. 
(IVGDR) was obtained in the 1937 



"?he first evidence of isovector giant dipole resonance 
l|. Extensive experimental and theoretical studies 
started about ten years later and provided simple expressions that systematically reproduce 
the excitation energy and the strength of the IVGDR in the whole nuclear chart [2i] . We are 
currently living in a sort of Renaissance era for this dipole mode since it plays a crucial role in 
the study of skin and halo structures in nuclei with large proton-neutron asymmetry, which 
can be obtained with radioactive nuclear beams. There is in fact experimental evidence 

1- known giant dipole resonance. 



The nature of this dipole 



that in nuclei with neutron excess, in addition to the we 
an accumulation of strength appears at low energies 
resonance is an open problem. In light nuclei, several theoretical calculations seem to relate 
it to nonresonant independent single-particle excitations of loosely bound neutrons [l2|, but 



there are some experimental results, for example for the ^^Ne [9'], which disagree with this 
picture. In medium and heavy nuclei this excitation seems to be collective, according to 
several theoretical calculations [12I, and it has been interpreted as a resonant oscillation of 
the neutron skin against the remaining isospin saturated neutron-proton core. It has been 
called pygmy dipole resonance (PDR). This denomination is often used to indicate this new 
collective mode, which in general exhausts a small fraction of the energy-weighted sum rule 
(EWSR). Sometimes, on the other hand, it is used in literature to generically indicate low- 
lying dipole excitations. The degree of collectivity and its evolution with the mass number 
is one of the most important open questions, but, as recently explained 13|], several other 



open problems remain in connection with this low-energy dipole excitation. One of them 
is the isoscalar or isovector nature. Theoretically, it can be investigated by calculating the 
neutron and proton transition densities; experimentally, it can be investigated by inducing 
the same excitation through different processes including electromagnetic processes, such 
as (7,7'), (e, e') or Coulomb excitation, on one hand, and strong processes, such as {a, a') 
or ions collision, on the other hand. Another important problem is the energy location of 
this excitation with respect to the neutron threshold, in connection with the question of 
whether some part of the PDR low-energy tail might be missing. The evolution of the PDR 
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with neutron-proton asymmetry as well as the existence of a connection between the PDR 
in stable and exotic nuclei are also of great interest. 

Here, we investigate these problems by focusing on neon isotopes and = 16 isotones 
through a fully consistent axially-symmetric-deformed quasiparticle random phase approx- 
imation (QRPA) calculation based on Hartree-Fock-Bogolyubov (HFB) states. The same 
Gogny DIS effective force |l4| has been used both in HFB and QRPA, which ensures the 
consistency of the calculations. This approach is essential in open-shell nuclei, where pairing 
correlations play an important role. It also allows the treatment of other isotopic and iso- 
tonic chains far from the closed major shells of stable nuclei. It has already been employed 
to study g,aut resonances of sUieon and magnes.um isotopes Q and very recently the heavy 
deformed ^ss-y jj^g|_ Exotic spherical ^'^Ni, ^"^^Sn, ^^^Sn, and ^''Spb were previously studied 
in the corresponding spherical HF+RPA calculation 17 |. 



Some stellar phenomena, such as r-process nucleosynthesis, are particularly sensitive to 
the low-energy tails of dipole responses [18|. Theoretical microscopic studies of photoab- 
sorption and radiative neutral capture cross sections pjj are essential in the regions of the 
nuclear chart where data are absent. All the models which pretend to be used in the whole 
nuclear chart without changing any parameters must be tested against the available exper- 



imental data. Also for this reason, after the calculations of 



15 



171 we focus here on neon 

n 

isotopes and = 16 isotones since experimental data on ^^Ne are available |9| and those on 



24/- 



O will appear soon [20|. Work on Ni is also in progress due to the experimental interest 



in the nature of the PDR 



2l| in this nucleus 10|. In connection to this, an interesting arti- 



cle 22[ relating the PDR to the nuclear symmetry energy 23| and to the isospin-dependent 
components of effective nuclear interactions recently appeared. 

Our paper is organized as follows: We present our model in Sec. II and our results in 
Sec. HI. Our work is summarized in Sec. IV. 



II. MODEL 

In this section, we briefly sketch our approach based on the QRPA on top of HFB calcu- 
lations. Details of the formalism, as well as of ^^~^'^Mg and ^^^^'^Si dipole responses can be 
found in [15]. The HFB equations are solved in a finite harmonic oscillator (HO) basis. As 
a consequence, the positive energy continuum is discretized. We have checked the stability 
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of the single-particle levels with the number of major HO shells. Finally, according to |15| . 
we considered a model space including 9 HO major shells. This is a large enough space 
for the neon isotopes and = 16 isotones considered here. To check it, we performed the 
calculation in a larger HO basis for ^^Ne. For 10 and 11 major shells, we obtain the same 
results as for 9 major shells both for spurious and low-lying states. All HFB quasiparticle 
states are used to generate the 2-quasiparticle (2-qp) excitations. This means that no cut 
in energy or in occupation probabilities is introduced. As already emphasized, we use the 
same nucleon-nucleon effective force, the Gogny DIS [14], both for HFB and QRPA calcu- 
lations in all particle- hole (ph), particle-particle (pp), and hole-hole (hh) channels. This is 
very important in order to avoid very dangerous inconsistencies. For this reason, since the 
Coulomb exchange field is not taken into account in HFB, the corresponding QRPA terms 
have been set to zero. 

According to the symmetries imposed in the present axially-symmetric-deformed HFB 
calculations in even-even nuclei, the projection K of the angular momentum on the symmetry 
axis and the parity tt are good quantum numbers. Consequently, QRPA calculations can be 
performed separately in each K"^ block. In an axially-symmetric-deformed nuclear system, 
the response function of a given J'^ contains different K"^ = 0^ ,±1'^ , ...,±J'^ components. 
In spherical nuclei, all these components are degenerated in energy, and then the response 
functions associated to any multi-polarity can be obtained from K'" = 0"^ results only. To 
solve the well-known (Q)RPA matrix equation 

I I , (1) 




where a;„ are the energies of the QRPA excited states, we use the same numerical procedure 



recently applied to study the giant resonances of the heavy deformed ^^^U 16|. It is based 
on a massive parallel master-slave algorithm. For a single solution of Eq. ([T]) the QRPA 
provides the set of amplitudes and Y"' describing the wave function of the excited state 
\n) in terms of the two quasiparticle excitations. Let us define A^2qp as the number of all the 
possible 2-qp excitations for a given K"^ block. The well-known normalization of the QRPA 
amplitude can then be written as 

E [(^2,p)' - iy2j'] = 1- (2) 

2qp=l 



From this equation, it is easy to isolate the relative neutron and proton contributions by 
summing separately over the two species. In the following, we indicate the corresponding 
percentage contribution with the notation %(z/) and %(7r). In order to measure the degree of 
collectivity of a specific excited state, we consider, inspired by [24], another index, called A^*, 



which represents the number of states with [(Xgqp)^ — (^2qp)^] — ^^^^ '^^^^ 



separate neutron (A^*) and proton (A^*) contributions. Note that in Ref. j2J] RPA, instead of 
QRPA, calculations are performed. In this case, A^2qp is replaced by A^ph, which represents 
the number of possible particle-hole excitations. The A^* value in the RPA case (called 
here N*^^^) represents the number of states with [{X^^)'^ — (l^p^)^] > In this case, 

when the excitation is produced by a single particle-hole state N*^^^ = 1 while in an ideal 
collective case all the particle-hole excitations contribute with the same statistical weight 
so that A^*^^^ = N*^^"^ + N^RPA = N^^. In the QRPA, the situation is a little different 
because 2-qp excitations imply the inclusion of hole-hole and particle-particle transitions. 
Anyway, as the A^* values become greater, the response becomes more collective. Obviously, 
the quantity A^ph, as well as A'2qp, is related to the size of the configuration space, which 
depends on the number of major oscillator shells one decides to consider and overall on the 
symmetry scheme (spherical, axially-symmetric-deformed,...) of the calculation. For this 
reason, a direct comparison of our results with the ones of Ref. 2J] is not possible. In Sec. 
nil Al we perform a comparison between this kind of analysis in RPA and QRPA schemes. 
Note that even if for spherical nuclei the results obtained for K = {] are strictly equivalent 
to the ones obtained for = 1, we choose to perform the analysis in the K = Q block. 
In this way, the number of the possible 2-qp configurations for all the nuclei we consider is 
always the same, precisely A"2qp=4832 (half obviously being neutron-neutron 2-qp excitations 
and half proton-proton 2-qp excitations). Because we always consider the same number of 
possible 2-qp excitations, it is not necessary to introduce other indexes, such as the ratio 

III. RESULTS 



The ground-state properties of the considered nuclei are obtained by solving HFB equa- 
iions. Axially symmetric HFB+DIS results for the whole nuclear chart can be found in 



25| . For deformation parameters, neutrons and protons root-mean-square radii, pairing 
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and Fermi energies we refer the reader to this web site. In Fig. [H we just plot the ground- 
state density profiles for the neon isotopes and = 16 isotones evaluated in this approach 
with 9 HO major shells. In the case of deformed nuclei, the drawn density is the spherical 
projected one, expressed as a function of the r spherical coordinate. For spherical nuclei the 
spherical radius r coincides with the radial distance (r_|_) from the symmetry axis z when the 
longitudinal position is = or, equivalently, r = when = 0. We observe in Fig. [1] 
the formation and increase of a neutron skin on the nuclear surface when the neutron-proton 
ratio increases. For ^^Ne a proton skin appears. 

Figures |2] and [3] show the B(E1) distributions for = 16 isotones and neon isotopes 
respectively. All the A^ = 16 isotones we consider are spherical, according also to prediction 



of A^ = 16 subshell closure from stability to the neutron drip line [26|]. Also, the analyzed 
neon isotopes, except ^^Ne (/3 = 0.5), 22Ne (/3 = 0.5) and ^^Ne (/3 = 0.2), turn out to be 
spherical. We write in each panel of Fig. |3]the axial deformation parameter /3 corresponding 
to the minimum of HFB potential energy. For ^^Ne the two minima at /3 = 0.2 and /3 = —0.2 
are quasidegenerated and the spherical (/3 = 0) configuration is just 100 keV less bound 
than the prolate one. For this reason, for ^^Ne both the responses (/3 = 0.2 and /3 = 0) 
are given in Fig. [31 In the spherical symmetry case, the K'"=^^ and |i^'|'^=l^ states 
are degenerate. In deformed ^°Ne, ^^Ne, and ^^Ne nuclei the strength splits up into two 
components corresponding to two different angular momentum projections K. As expected, 
since these nuclei are prolate, the GDR peak at low energy ~ 20 MeV) corresponds to 
the K=^ components. The l-ft'Hl GDR strength is concentrated around 25 MeV. For ^°Ne 
and ^^Ne nuclei, the lOn < 14 MeV excitations are negligible. In all the other cases one or 
some peaks appear in this low-energy region. The giant resonance is located around 20-25 
MeV. 

Each spectrum has been cleaned by removing the corresponding translation spurious 
state. In order to identify it, a small renormalization factor a of the residual interaction has 
been introduced. Figure H] shows the evolution of the first ten Un eigenvalues as a function of 
this a parameter. As clearly shown, only the spurious state is affected by this factor while 
all the other states, including the ones corresponding to the low-energy peaks, stay quite 
constant. We checked also the stability of the B(E1) values with a in the case of ^^Ne, and 
the strengths of all states but the spurious one remain always the same. The spurious state, 
which turns out to be easily identifiable, is thus decoupled from the physical spectrum. 
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A. ^^Ne and = 16 isotones 

We start our analysis from ^^Ne, which belongs both to Z = 10 and = 16 chains. The 
low-energy dipole excitations of this nucleus have been experimentally studied jol- Several 



theoretical calculations on this nucleus have been already done 27H3l|. In Table [T] we 
present a detailed analysis of the microscopic structure of the low-lying dipole modes in 
^^Ne (and in = 16 isotones) considering the first low-energy state. Instead to show the 
contribution of the main two quasiparticle excitations, we prefer to give the results in terms 
of the corresponding particle-hole (and particle-particle) transitions identified by the usual 
spectroscopic notation. The energies corresponding to each single-particle transition for each 
nucleus are the values given between parentheses in Table [B Focusing at first on ^^Ne, it 
appears that the state at ujn = 10.7 MeV is dominated by the two quasiparticle excitations 
corresponding to the u 2si/2— ^2p3/2 transition. In terms of percentages this contribution 
turns out to be of 67.6 %, a value very close to the one of {sol. According to [2^, the second 
contribution (9.5 %) arises from the u 2si/2—^2pi/2 excitation. Two-quasiparticle proton 
excitations are 8.1 %. An important feature of our results is an appreciable contribution 



(8.7 %) of the u ld5/2— ;'lf7/2 configuration. The experimental results of {q] seem to suggest 
that the low- lying dipole excitations around 6-10 MeV involve more transitions with respect 
to those only characterized by 2s^y^2 states. As stressed in [qI theoretical private calculations 
in the approach of [28] predict a nearly equal weight of z/ 2si/2— ;'2pi/2 and z/ ld5/2— !'lf7/2 
transitions in the low lying = 1~, which turns out to be dominant in this model. In 
connection with the experimental results of [91], also the structure of the second main peak, 
at Un = 17.9 MeV is interesting. It is mainly generated by the u lpi/2— ^ld3/2 (49.3 %), u 
ld5/2^1f7/2 (24.4 %) and u ld5/2->2p3/2 (14.42 MeV) (13.1 %) transitions. 

Concerning the A^ = 16 isotones, the other nuclei we consider here are ^^O, ^^Mg, and 
^"Si. While the low-energy peak is expected in ^^O, one can observe in Fig. [2] that an 
equivalent structure appears at Un — 12 MeV also for ^®Mg and ^"Si in spite of their small 
neutron-proton asymmetry. As expected and as shown in Fig. [5] the centroid energy of the 
low- lying El states increases with the decrease of neutron excess and the B(E1) value of the 
first peak is suppressed. On the other hand, the total strength up to 14 MeV decreases very 
slowly from ^^Ne to ^°Si. 

Turning to the microscopic analysis of Table [H one can observe that the structure of the 
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first peak in ^ O is quite similar to the one of ^^Ne except for the absence in ^ O of the v 
ld3/2— )-2p3/2 contribution, as expected for a particle-particle configuration in a closed-shell 
nucleus. This particle-particle transition turns out to be very important in ^^Mg, being of 
the same order of magnitude as the main particle-hole transition v 2si/2— '■2p3/2- Finally, in 
^"Si the contribution of proton excitations becomes appreciable, being of the order of 30%. 
In the three other cases, the neutron contribution to the low energy peak is more than 90%. 
Even if a detailed study of the GDR region is beyond the aim of this paper we consider here 
and in the following the behavior of the main GDR peak in order to compare it with the one 
of low-lying states. For this reason in Table [T] the total proton and neutron contributions 
to the main GDR peak are reported. The difference with respect to the corresponding first 
peak clearly appears. It is interesting to observe that even in the GDR peak the neutron 
contribution is larger than the expected N j A value in the cases of ^^O and ^^Mg nuclei. 
Note always from Table [T] that if instead of treating only the main GDR peak we consider 
the total proton and neutron contribution arising from all the states with 15 MeV < < 
30 MeV and B(E1) > 0.5 e^ fm^, the result is quite similar. 

Further remarks about the collective behavior of the analyzed dipole excitations in these 
= 16 isotones arise from the calculation of the A^*, A*"* and A^* indexes defined in Sec. 
inland reported in Table HTl together with the isospin percent contributions %{y) and %(7r). 
Because we start the analysis from ^'^O, some comments are needed. The little values of A^* 
numbers (even in the GDR case) with respect to A"2qp must not be surprising. This apparent 
disproportion is related to the huge number A'2qp = 4832, which is due to the choice of 
performing calculation in QRPA, including in this way not only the particle-hole transitions 
but also the particle-particle and the hole-hole ones without any cut in excitation energy. 
It is important to observe that in the corresponding spherical RPA calculation (which for 
^^O is strictly equivalent to the QRPA one) the number of possible particle-hole excitation 
is largely reduced (A'2qp = 4832 — i- A^ph = 64), increasing the relative importance of A^* 
indexes, which in this case become A^*-^^^=3 and N*^^^=2 at a;n=9.1 MeV and N*^^^=7 
and N*^^'^=3 at u;„=20.5 MeV. 

Coming back to the analysis of the QRPA results shown in Table [Tll one can observe that 
the number of state which contribute with a weight greater than as defined in Sec. 

IIT| lightly increases from ^^O to ^^Si both for low-energy mode and for the GDR. The A"* 
index generally decreases with Z while A^* increases. The number of configurations which 
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contribute to the low energy mode is about the 60% of the configurations involved in the 
GDR. This value reduces to about 40% just in the ^^Mg case. This apparently less collective 
behavior of the 11.6 MeV peak is compensated by the appearance of another quite close 
peak at 11.8 MeV, which is dominated by the v ld3/2— ;"2p3/2 (62%) and the v 2si/2— ^2p3/2 
(24%) transitions. 

In Fig. [6] we display the neutron and proton transition densities bp for ^^Ne and for 
the other = 16 isotones for the low-energy and the main GDR peaks. Note that since 
the nuclei we considered turn out to be spherical, even if we performed our calculation in 
cylindrical coordinates, we can and prefer to present the transition densities as a function of 
r instead that as a function of and r^. The lower panels of FigE] refer to the main GDR 
peak. They display for all the = 16 isotones a radial dependence that is characteristic for 
the isovector GDR: The proton and neutron densities oscillate with opposite phases. The 
low-energy case, shown in the upper panel, is more complex but is quite similar in all the 

= 16 isotones. It presents an isovector character in the center of the nucleus (for r < 2 
fm) and an isoscalar one in the surface region (r ~ 3 fm). For r > 4 fm one can observe 
an oscillating finite contribution of the neutrons and a negligible one of the protons. The 
spatially extended structure of the 2si/2 state is responsible for the shape of the neutron 
transition density. In order to illustrate this point we show in Fig. [7] the ^^O neutron 
transition density for the state at 9.1 MeV, as already plotted in the first panel of Fig. |6l 
and the unperturbed transition densities for v 2si/2— ;'2p3/2 (8.92 MeV) and v 2si/2— ^►2pi/2 
(10.33 MeV) excitations. It clearly appears that the behavior at the center of the nucleus is 
quite similar, refiecting the crucial role of the v 2si/2 state, while in the surface region the 
two unperturbed transition densities differ from the one obtained in QRPA, reflecting the 
important role of correlations in this region. For r > 5 fm, the three transition densities 
turn to be similar. 

The result obtained for the transition densities of the low-energy excitations of A^ = 16 
isotones does not coincide with the usual isoscalar rep resentation of the PDR obtained in 



many theoretical calculations but for heavier nuclei 12| . Transition densities of the low- lying 
dipole states in ^^Ne can be found in 30 1 and in j28|, while results for ^"^O are presented in 



32| and 2J| where a phenomenological RPA is used. In 32|, the isoscalar character of the 



excitation inside the nucleus clearly appears while the isovector behavior of the analyzed 
low energy excitation (however not the lowest one) in 2M, is interpreted as the tail part of 
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the GDR. At the moment, the results obtained for the transition densities of the low- lying 
dipole state in these nuclei seem to be model dependent, and the problem remains open. 

To summarize, the low-lying dipole state of = 16 isotones possesses a small but finite 
collective nature. While the strength of this state increases with the neutron-to-proton 
ratio, the degree of collectivity stays quite constant. These considerations, as well as the 
study of the transition densities, suggest that some interesting low-energy properties can be 
experimentally obtained also from the less exotic ^^Mg and ^°Si nuclei. 



B. Neon isotopes 

The analysis of neon isotopes proceeds along the same lines. Figure |8] is the equivalent 
of Fig. [5] but for the Z = 10 instead of = 16 chain. Also in this case, the energy of the 
low- lying peak decreases with the isospin asymmetry and the B(E1) strength increases. As 
shown in the lower panel of FigJHlby comparing the B(E1) value of the low-energy peak with 
the value obtained summing the B(E1) up to 14 MeV the fragmentation increases with the 
number of neutrons. It merges also from Fig. [3l 

The lower mass isotope, which presents a non-negligible strength, is ^^Ne, as recently 
found also in 3l|. According to the HFB ground-state calculation [25], this nucleus turns out 
to be lightly deformed. Nevertheless, owing to the soft deformation with prolate (/3 = 0.2) 
and oblate (/3 = —0.2) configurations mostly degenerate, we prefer to consider ^"^Ne as 
spherical in our microscopic analysis. The comparison with other neon isotopes becomes in 
this way more easy and more consistent. Anyway, in Fig. [3]we reported the B(E1) results for 
^^Ne both in /3 = 0.2 and (3 = cases. The deformed response is obviously more fragmented 
with the splitting of A' = and A' = 1 components but the correspondences of low-energy 
and GDR excitations in /3 = 0.2 and /3 = cases are evident. For example, focusing on 
the little low-energy excitation, one can observe that the single peak at co'„=13.6 MeV with 
B(E1)=4.5 10~^ e^ fm^ in the spherical case is split, in the prolate case, in one peak for the 
A' = component at a;„=13.1 MeV with B(E1)=0.8 10"^ e^ fm^ and in three A' = 1 peaks 
at w„=13.4 MeV, w„=13.8 MeV and w„=13.9 MeV with B(E1) = 1.6 10^^ g2 ^^^2^ B(E1)=0.7 
10~^ e^ fm^ and B(E1)=2.1 10~^ e^ fm^, respectively. Turning to the microscopic analysis in 
the spherical case, the peak at 13.6 MeV is mostly due (65.4%), as in ^^Ne, to u 2si/2— ^2p3/2 
(15.03 MeV) transition. In this case, the u 2si/2— !'2pi/2 (15.93 MeV) transition turns out to 
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be more relevant (25.7%) with respect to ^^Ne. Two other transitions, v ld3/2— )'2pi/2 (19.78 
MeV) and tt lpi/2— ;"2si/2 (16.14 MeV), are of the order of 2%. Each of the remaining ones 
contributes individually less than 1%. The transition densities of this state, shown in Fig. 
O are quite interesting revealing a neutron fluctuation larger than the proton one. The very 
small proton fluctuation is anyway in phase with the neutron one, reflecting the isoscalar 
nature of this resonance. In contrast, the transition densities corresponding to all the other 
remarkable peaks shown an evident isovector behavior typical of the GDR. Furthermore, as 
deduced from A'^^ and A^* indexes reported in Table UTTt the excitation at 13.6 MeV involves 
principally neutron configurations and is little more collective than the successive one at 
16.4 MeV. 

Since we have already discussed in the previous subsection ^^Ne, we turn now to ^®Ne. 
In this case, the dipole strength distribution is more fragmented. The first two peaks 
are essentially of single-particle nature. They are related to the transitions from the less 
bounded ld3/2 neutron state: v ld3/2^2p3/2 (8.26 MeV) (90.4%) for the 8.1 MeV peak and 
V ld3/2^2pi/2 (10.00 MeV) (83.1%) for the second one at 9.5 MeV. The corresponding A^* 
and A^* indexes shown in Table IIIII are lower with respect to the ones related to higher 
energy excitations. The state at 11.1 MeV corresponds to the low energy peak of ^^Ne. The 
dominant component is the same, that is, v 2si/2— >2p3/2 (10.82 MeV) (67.2%). The degree 
of collectivity is very similar, according to Table llllt as well as the neutron and proton tran- 
sition densities illustrated in Fig. [TOl The successive peak at 13.7 MeV is really interesting. 
As reported in Table IIIH the number of 2-qp neutron excitations related to this state is 
appreciably greater than the one of the lower states and even of the successive high-energy 
state. In this 13.7 MeV peak, the transition from ld3/2 neutron state turns to play a ma- 
jor role: 68.8 % from v ld3/2— >lf5/2 (14.09 MeV). The prominent role of transitions from 
the weakly bound ld3/2 neutron state is reflected in neutron densities, no more peaked at 
the interior of the nucleus as for the transitions involving the v 2si/2 state. The isoscalar 
behavior of the state dominated by these v ld3/2 transitions is also evident. 

Results obtained for ^°Ne are very similar to those for ^'^Ne concerning both the micro- 
scopical analysis and the transition densities shown in Fig. [TT] The first two peaks are 
still less collective: 93.5% from v ld3/2— '■2p3/2 (6.57 MeV) in the first peak at 6.6 MeV 
and 89.9% from v ld3/2— )'2pi/2 (8.09 MeV) in the second one at 7.9 MeV. The respective 
A"* and A^* indexes are lower than the corresponding ones in ^^Ne. The third peak at 10.5 
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MeV corresponds again to the low-energy peak of ^^Ne with again a similar microscopical 
structure: it is dominated (67.4%) by the v 2si/2— )-2p3/2 (10.20 MeV) transition. This peak 
is characterized by the same 2-qp configurations and the same transition densities as the 
= 16 isotones and the neon isotopes ^^Ne and ^'^Ne. According to the value of A^* it 
seems to be little more collective with respect to the corresponding excitation in the other 
nuclei. In the fourth peak at 12.8 MeV the main transition v ld3/2— ;'lf5/2 (12.08 MeV) (54.3 
%) is more mixed to the v ld5/2^1f7/2 (12.87 MeV) (18.6 %) one. As in ^SNe, this fourth 
peak is again the more collective state with the typical isoscalar behavior of the transition 
densities. The features of this fourth peak of ^®Ne and ^°Ne are more in touch with the 
standard representation of the PDR. 

C. Proton electric dipole resonance in ^^Ne 

As can be observed in Fig. [1], the ^®Ne exhibits a proton skin in spite of the presence 
of the Coulomb barrier. It is interesting to investigate if in this case a low-energy proton 
electric dipole resonance appears, paying attention to its degree of collectivity. Until now, 
only a few studies of dipole excitations in proton-rich nuclei have been done. In a work 



based on relativistic QRPA on A^ = 20 isotopes and Z = 18 isotones [33|, it has been shown 
that proton pygmy dipole resonances may appear when approaching the proton drip line. 
Here we focus on the ^^Ne case. In this nucleus the proton Fermi level is ld5/2. In the 
conditions of our HFB ground-state calculations (9 HO major shells, absence of Coulomb 
exchange field) the occupation probability of this level is i;^=0.32 while the corresponding 
single particle energy is CTridg/j = —2.14 MeV. The root-mean-square radius calculated for 



this orbital is a /(?"^)7rid5/2 = 3.40 fm, to be compared to the total proton root-mean-square 



radius a/ (r^)^ = 2.84 fm. 

Turning to the excited states, the first panel of Fig. [3] shows the ^^Ne B(E1) distribution 
which is dominated by the GDR around 23 MeV. On the other hand, two distinct peaks 
appear below 15 MeV. Their microscopical structures are given in Table IIVI The first 
peak at 14.2 MeV is essentially due to proton excitations. The dominant configurations 
correspond to transitions from the weakly bound proton state ld5/2. Contributions of 10% 
order of magnitude arise from transitions from lpi/2, 2si/2 and lp3/2 proton states. The 
total neutron contribution is less than 8%. The situation is different for the second low- 
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energy peak, located at 14.8 MeV, which is mainly due to neutron and proton lpi/2— ?• 2si/2 
transitions. The behavior of the first peak, characterized by a superposition of many 2-qp 
proton configurations, is specular with respect to one of the other neon isotopes and = 16 
isotones: Only the roles of neutron and proton are obviously inverted. This clearly appears 
also by looking at the transition densities of Fig. [T21 In this case, the proton and neutron 
are in opposite phase in the nuclear interior and in phase in the surface region, where the 
proton contribution is larger than the neutron one. On the other hand, the 14.8 MeV excited 
state already displays a pronounced isovector behavior as the GDR here appearing around 
23 MeV. It is also very important to observe that the degree of collectivity of the first low 
energy state is the same as the one of the GDR since the A^* values (shown in Table IIIip 
are the same. The A^* value is greater for the first low-energy peak revealing once again the 
dominant role of protons in this state. The correspondence between the results obtained 
for the low-lying excitation of this proton-rich nucleus and the neutron-rich nuclei here 
analyzed suggests that systematic studies of proton-rich nuclei in this mass region should 
be performed. 

IV. SUMMARY AND CONCLUSIONS 

We have studied the dipole excitations of neon isotopes and = 16 isotones focusing 
in particular on the low-excitation-energy region. We performed a fully consistent axially- 
symmetric-deformed HFB+QRPA calculation using the DIS Gogny nucleon-nucleon effec- 
tive force. For the A^ = 16 isotones we considered, we obtained excitation modes between 
9 and 12 MeV. All these modes presented a small but finite collective behavior: also other 
transitions with respect to the dominant v 2si/2— '■2p3/2 contribute to generate the low-lying 
resonance. On the other hand, the spatial structure of the less-bound neutron state, that 
is, the V 2si/2, is the main responsible of the neutron transition density and this is reflected 
in the nontrivial macroscopic isoscalar or isovector behavior of the low energy resonance. It 
holds also for the neon isotopes where the behavior of the transition densities is essentially 
due to spatial structure of the main hole state contributing to the particle-hole transitions. 
Furthermore, also for the neon isotopes, we found that the strength of the low-energy state 
increases with the neutron-to-proton ratio. By increasing the number of neutrons the dipole 
strength distribution becomes more fragmented. In ^^Ne and ^*^Ne the occupation of the v 
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lds/2 shell leads to the appearance of another state just above the one corresponding to the 
excitation belonging to all the = 16 nuclei analyzed. This state is more collective and it is 
characterized by an isoscalar behavior of the transition densities, revealing features that are 
more in touch with the standard representation of the PDR characterizing heavier nuclei. 

Finally, for the proton-rich ^^Ne, we obtained a low-energy resonance at 14.2 MeV whose 
transition density behavior is quite similar to the one of the low- lying excitation of N — 16 
isotones once the role of neutrons and protons is inverted. The degree of collectivity of this 
state is the same as the one of the main GDR peak of this nucleus. 
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240 


26Ne 


28Mg 


30Si 


First peak 


a;„=9.1 McV 


C^n = 10.7 McV 


tj„=11.6 MeV 


tj„=12.2 MeV 


V 2si /o— >2pQ/o 


73.5 % (8.92 MeV) 


67.6 % (10.52 MeV) 


39.7 % (11.68 MeV) 


43.8 % (12.61 MeV) 


V ld3/2^2p3/2 


0.0 % (9.26 MeV) 


2.8 % (10.82 MeV) 


32.7 % (11.93 MeV) 


7.0 % (12.79 MeV) 


V 2si/2— >-2pi/2 


10.0 % (10.33 MeV) 


9.5 % (12.44 MeV) 


5.8 % (13.98 MeV) 


6.6 % (15.23 MeV) 


V ldc;/o — ^-Vii 1^ 


8.7 % (12.87 MeV) 


8.7 % (13.68 MeV) 


8.2 % (14.18 MeV) 


11.1 % (14.55 MeV) 


V Ids/o— 7'2pi 11 


0.0 % (10.67 MeV) 


0.1 % (12.74 MeV) 


2.7 % (14.22 MeV) 


1.1 % (15.41 MeV) 


V Ipi /O^ldQ/O 


1.2 % (17.70 MeV) 


1.1 % (18.50 MeV) 


0.5 % (19.10 MeV) 


0.3 % (19.60 MeV) 


V total contribution 


94.6 % 


91.9 % 


90.6 % 


70.9 % 


TT lpi/2~^2S]^/2 


2.0 % (15.19 MeV) 


4.1 % (15.97 MeV) 


5.6 % (16.22 MeV) 


18.5 % (15.92 MeV) 




2.7 % (15.89 MeV) 


2.0 % (17.60 MeV) 


0.8 % (18.97 MeV) 


0.9 % (17.15 MeV) 




(lp3/o-> Idc/o) 


(lp3/o— ldc;/n) 


(lp3/9^ ldc;/9) 


(ldr,/9^2po/r,) 


TT ldc;/9— >-lf7/9 
Of z ' / ^ 


0.0 % (24.80 MeV) 


0.9 % (17.48 MeV) 


1.9 % (16.02 MeV) 


7.7 % (14.33 MeV) 


TT total contribution 


5.4 % 


8.1 % 


9.4 % 


29.1 % 


Main peak (GDR) 


a;„=20.5 MeV 


cjn=21.9 MeV 


a;n=22.6 MeV 


a;n=23.1 MeV 


u total contribution 


80.5 % 


64.5 % 


78.5 % 


54.4 % 


TT total contribution 


19.5 % 


35.5 % 


21.5 % 


45.6 % 


Peaks with 


B(E1) > 0.5 e^ fm^ 


and 15 MeV 


< < 


30 MeV 


u total contribution 


81.6 % 


78.0 % 


82.9 % 


52.7 % 


TT total contribution 


18.4 % 


22.0 % 


17.1 % 


47.3 % 


N/A 


66.7 % 


61.5 % 


57.1 % 


53.3 % 



TABLE I: Main particle-hole (and particle-particle) configurations contributing to the first low- 
energy dipole excitations of AT = 16 isotones. Total neutron and proton contributions are also 
given for the low-energy peak, for the main one corresponding to the GDR, and from the sum of 
all the states with 15 MeV < a;„ < 30 MeV and B(E1) > 0.5 e^ fm^. 
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ojn (MeV) A^* 


m 


m 

TV 


%M % (vr) 












9.1 


44 


34 


10 


94.6 % 5.4 % 


20.5 


74 


62 


12 


80.5 % 19.5 % 


26Ne 










10.7 


48 


30 


18 


91.9 % 8.1 % 


21.9 


82 


56 


26 


64.5 % 35.5 % 


28Mg 










11.6 


46 


24 


22 


90.6 % 9.4 % 


22.6 


114 


82 


32 


78.5 % 21.5 % 


30Si 










12.2 


54 


32 


22 


70.9 % 29.1 % 


23.1 


88 


54 


34 


54.4 % 45.6 % 



TABLE II: Values of the collectivity indexes for the low-lying dipole excitation and GDR in the 
N = 16 isotones. 
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, (MeV) A^* 


N* 


N* 




%(7r) 


i^Ne 












14.2 


72 


12 


60 


7.8 % 


92.2 % 


14.8 


46 


10 


36 


44.3 % 


55.7 % 


23.0 


72 


16 


56 


21.9 % 


78.1 % 


24Ne 












13.6 


52 


44 


8 


97.7 % 


2.3 % 


16.4 


48 


32 


16 


55.1 % 


44.9 % 


23.9 


84 


60 


24 


65.0 % 


35.0 % 


26Ne 












10.7 


48 


30 


18 


91.9 % 


8.1 % 


21.9 


82 


56 


26 


65.5 % 


35.5 % 


28Ne 












8.1 


36 


24 


12 


98.7 % 


1.3 % 


9.5 


40 


26 


14 


98.4 % 


1.6 % 


11.1 


46 


30 


16 


95.2 % 


4.8 % 


13.7 


86 


72 


14 


97.1 % 


2.9 % 


16.1 


78 


32 


46 


5.0 % 


95.0 % 


18.3 


102 


76 


26 


91.5 % 


8.5 % 


20.9 


124 


92 


32 


70.2 % 


29.8 % 


30Ne 












6.6 


28 


18 


10 


98.9 % 


1.1 % 


7.9 


34 


20 


14 


98.9 % 


1.1 % 


10.5 


56 


40 


lb 


96.9 70 


10/ 
O.l 70 


12.8 


80 


64 


16 


97.0 % 


3.0 % 


16.0 


68 


26 


42 


3.2 % 


96.8 % 


18.7 


142 


98 


44 


68.3 % 


31.7 % 



TABLE III: Values of the collectivity indexes for the main excitations in the neon isotopes. 
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First low-energy peak 


Second low-energy peak 


Main peak (GDR) 




tj„=14.2 MeV 


a;„=14.8 MeV 


Wn=23.0 MeV 


TT ld5/2^1f7/2 (15.99 MeV) 


29.1 % 


1.6 % 




TT ld5/2^2p3/2 (15.98 MeV) 


23.2 % 


6.3 % 




TT lpi/2^2si/2 (13.75 MeV) 


11.7 % 


35.2 % 




TT 2si/2^2p3/2 (18.19 MeV) 


9.3 % 


1.2% 




TT lp3/2^1d5/2 (16.95 MeV) 


8.5 % 


9.3 % 




TT 2si/2-)-2p3/2 (19.14 MeV) 


2.7% 


0.5% 




TT total contribution 


92.2 % 


55.7 % 


78.1 % 


V lp3/2^1d3/2 (15.78 MeV) 


3.8 % 


2.7% 




V lpi/2^ 2si/2 (14.18 MeV) 


3.6 % 


41.0 % 




V total contribution 


7.8 % 


44.3 % 


21.9 % 



TABLE IV: Main particle-hole (and particle-particle) configurations contributing to the first two 
low-energy peaks of ^^Ne. Total neutron and proton contributions are also given both for the 
low-energy peaks and for the one corresponding to the GDR. 
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FIG. 1: (Color online) Neutron (solid lines) and proton (dashed lines) ground-state density profiles 
for neon isotopes (left panel) and = 16 isotones (right panel) obtained through axially-deformed 
HFB calculation with DIS Gogny interaction. 
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FIG. 2: (Color online) B(E1) distributions for = 16 isotones in logarithmic scale. 
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FIG. 3: (Color online) B(E1) distributions for neon isotopes in logarithmic scale. For deformed 
isotopes ■^''■^^'■^'^Ne, the K = and K = ±1 components are drawn in black (blue online) and gray 
(green online) respectively. In the spherical (/3 = 0) case this splitting does not appear, and hence 
the K = and K = ±1 coincide (red online). 




FIG. 4: (Color online) Evolution of the first ten ujn eigenvalues with the a parameter for the neon 
isotopes and = 16 isotones, which present low-lying dipole excitations. The lowest energy state 
is the spurious one. It is plotted with continuous line (red online) as well as the state corresponding 
to the analyzed excitation (blue online). The remaining states are plotted with dashed black lines. 
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FIG. 5: (Color online) Energy of the low-lying peak (upper panel) and B(E1) transition strengths 
(lower panel) for the = 16 isotones as functions of mass number. 
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FIG. 6: (Color online) Neutron (full lines) and proton (dashed lines) transition densities for the 
24q^ -^^Ne, ^^Mg and ^'^Si nuclei. The corresponding excitation energies are given in the panels. 
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FIG. 7: (Color online) ^'^O neutron transition density for the first low-energy peak at 9.1 MeV 
calculated in QRPA (full line) compared to the unperturbed transition densities of v 2si/2— >'2p3/2 
(dot-dot-dashed line) and u 2si/2— ^2pi/2 (dot-dash-dashed line) excitations. 
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FIG. 8: (Color online) Energy of the first low-energy peak with B(E1)> 0.01 e^ fm^ (upper panel) 
and B(E1) transition strengths (lower panel) for the neon isotopes as functions of mass number. 
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FIG. 9: (Color online) Neutron (full lines) and proton (dashed lines) transition densities for ^^Ne. 
The corresponding excitation energies are given in the panels. 
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FIG. 10: (Color online) Neutron (full lines) and proton (dashed lines) transition densities for ^^Ne. 
The corresponding excitation energies are given in the panels. 
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FIG. 11: (Color online) Neutron (full lines) and proton (dashed lines) transition densities for '^'^Ne. 
The corresponding excitation energies are given in the panels. 
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FIG. 12: (Color online) Neutron (full lines) and proton (dashed lines) transition densities for ^^Ne. 
The corresponding excitation energies are given in the panels. 
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